LOWER ORDER TERMS FOR THE 
ONE-LEVEL DENSITY OF ELLIPTIC 
CURVE L-FUNCTIONS 



D. K. Huynh, J. P. Keating and N. C. Snaith 
School of Mathematics, 
University of Bristol, 
Bristol BS8 ITW, UK 

November 14, 2008 



Abstract 

It is believed that, in the hmit as the conductor tends to infinity, correlations 
between the zeros of elliptic curve _L-functions averaged within families follow 
the distribution laws of the eigenvalues of random matrices drawn from the 
orthogonal group. For test functions with restricted support, this is known 
to be the true for the one- and two-level densities of zeros within the families 
studied to date. However, for finite conductor Miller's experimental data reveal 
an interesting discrepancy from these limiting results. Here we use the L- 
functions ratios conjectures to calculate the 1-levcl density for the family of even 
quadratic twists of an elliptic curve L- function for large but finite conductor. 
This gives a formula for the leading and lower order terms up to an error term 
that is conjectured to be significantly smaller. The lower order terms explain 
many of the features of the zero statistics for relatively small conductor and 
model the very slow convergence to the infinite conductor limit. However, our 
main observation is that they do not capture the behaviour of zeros in the 
important region very close to the critical point and so do not explain Miller's 
discrepancy. This therefore implies that a more accurate model for statistics 
near to this point needs to be developed. 



1 Introduction 

The conjecture that the limiting statistical properties of the zeros of L-functions 
may be modeled by those of the eigenvalues of random matrices goes back to Mont- 
gomery |Mon73] ■ who introduced it in the context of the Riemann zeta-function. 
For the Riemann zeros this conjecture is supported by extensive numerical |Odl97| 
and theoretical |Mon731 |Hej94t IBK951 IBK96bl [RS96] calculations. The generaliza- 
tion to zero statistics within families of L-functions was developed by Katz and 



Sarnak [KS99al IKS99bj . and again there is much evidence supporting it [RubOlj . 
Random matrix models for the moments of the Riemann zeta-function on its crit- 
ical line and for central values of L-functions within families were introduced by 
Keating and Snaith [KSOObl IKSOOaj . and have since been developed extensively 
[CFOnl ICFK+nHl \GBKQ7[ IBJOTI IBJOHI ICFK+n8| . For more background, see [HiQS]. 

The random-matrix moment conjectures extend naturally to ratios of L-functions. 
The L-functions ratios conjectures were stimulated by the work of Farmer, who, in 
1995, made a conjecture for shifted moments of the Riemann zeta-function [Far95] . 
Nonnenmacher and Zirnbauer [NZj found formulas for the ratios of characteristic 
polynomials of random matrices coming from one of the classical compact groups. 
This was formalised and written up by Conrey, Farmer and Zirnbauer [CFZb] and 
lead to the development of corresponding ratios conjectures for L-functions in num- 
ber theory |CFZaj . 

The Birch/Swinnerton-Dyer conjecture asserts that the rank of an elliptic curve 
is equal to the order of vanishing at the central point of the associated L-function. 
The idea of using random matrix theory to predict the frequency of non-zero rank 
in families of elliptic curves was introduced by Conrey, Keating, Rubinstein and 
Snaith [CKRS021 ICKRS06] . An interesting extension of this is to find a random 
matrix model for elliptic curve L-functions of a given order of vanishing at the 
critical point. The first steps in this direction have been taken by Snaith |Sna05j 
and Miller/Dueiiez [Mil06| . but it is clear from Miller's numerical computations that 
there is a still simpler problem concerning the zero statistics of families of rank zero 
curves that is far from being understood. This problem is the main motivation for 
the work we shall report on here. 

According to the Katz/Sarnak philosophy [KS99at IKS99bj . zeros of families of 
L-functions show the same statistical behaviour as eigenvalues of random matrices 
drawn from one of the classical compact groups. The zeros of a family of elliptic 
curve L-functions with even (odd) functional equation should follow the distribution 
laws of eigenvalues of the even (odd) orthogonal group. Rigorous calculations |Mil021 
IMil04t lYouOG] show that as the conductor (the parameter that orders L-functions 
within a family) tends to infinity, the one- and two-level densities do indeed tend to 
the expected orthogonal forms for several different families of elliptic curves. That 
is, as the conductor tends to infinity, the zero statistics approach the scaling limit 
for large matrix size of the corresponding statistic for the eigenvalues of matrices 
from SO{2N) or SO{2N + 1). (Similar agreement with random matrix theory is 
shown for many other families of L-functions, see for example [DM06t IFI03t IGiilOSl 
IHR031 IHM071 IlLSOOl l6S99l iRRbl |Roy011 IRubOlj .1 The test functions involved in 
these calculations have a limited range of support, but nonetheless the evidence is 
compelling. Thus it was surprising to see in Miller's numerical results |Mil06| a 
distinct repulsion of the zeros from the central point for a family of L-functions 
of rank elliptic curves, because no repulsion is seen in the statistics of S0{2N) 
eigenvalues. Of course, in numerical computations the conductor is finite, and so it 
is clear that an explanation is needed for finite conductor statistics and how they 
approach the limiting S0{2N) statistic. 
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We do have a relatively complete understanding of the way in which the random 
matrix limit is approached for the zero statistics of the Riemann zeta function at 
a height T up the critical line as T ^ oo. Berry first wrote down an approximate 
formula describing the finite-T corrections to the random matrix limiting form for a 
statistic related to the 2-point correlation function in [Ber88j and showed that this 
described Odlyzko's data remarkably accurately. Later, a formula that is believed to 
capture all of the essential features was derived by Bogomolny and Keating [BK96aj . 
The terms in the Bogomolny-Keating formula that describe the corrections to the 
random matrix limit are often referred to as lower order terms. See [BK99] for 
an overview and numerical illustrations. More recently, Conrey and Snaith [CS071 
ICSj have shown how the Bogomolny-Keating formula and its extension to all n- 
point correlation functions can be recovered from the L-functions ratios conjectures 
[CFZaj . There have also been investigations of lower order terms in the zero statistics 
of various families of L-functions [FlOSl iMiIbl iMiIal [RRal lYouOSj . In particular, 
Conrey and Snaith have shown how such terms can also be recovered from the 
ratios conjectures |CSC)7] . It is thus natural in this context to seek the explanation 
for the surprising discrepancy observed by Miller in these lower order terms. 

In this paper we examine lower order terms in the 1-level density of the zeros 
of a family of elliptic curve L- functions. Specifically, we investigate even quadratic 
twists of an elliptic curve L-function, for which we calculate the zeros numerically 
with Rubinstein's Icalc [Rub]. Using the ratios conjectures we derive a formula for 
the 1-level density that describes convincingly the intricate structure of the numer- 
ical data away from the central point and so explains the rate of approach to the 
random matrix limit in this region. However, most interestingly, our formula fails 
to describe the region very close to the central point. To illustrate our main results, 
we plot in figure [1] a numerical evaluation of the 1-level density together with our 
formula. Miller's discrepancy corresponds to the region near to the origin. Our main 
conclusion here is then that the explanation for the zero distribution in this region 
lies beyond the models combining random matrix theory and arithmetical lower 
order terms considered so far; that is, these formulae are not sufficient to explain 
the discrepancy. We plan to explore augmented models that build on the present 
calculation to explain the phenomenon in a future paper with E. Duehez and S. J. 
Miller 

2 The 1-level density formula 

Let the L-function Le{s) associated with an elliptic curve E be given by the Dirichlet 
series 



where the coefficients (A(n) = a{n)/ ^/n, with ap = p+1 — #£'(Fp), #E(¥p) being the 
number of points on E counted over Fp) have been normalised so that the functional 




(2.1) 



n=l 
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Figure 1: 1-level density of unsealed zeros from up to height 0.6 of even quadratic 
twists of Le^i with < d < 100,000 for left and < d < 400,000 for right hand 
side, prediction (sohd), from ()2.52p . versus numerical data (bar chart). 



equation relates s to 1 — s: 

Here M is the conductor of the elliptic curve E; we will consider only prime M. 
Also, uj(E) is +1 or —1 resulting, respectively, in an even or odd functional equation 
for Le- 

Let LE{s,Xd) denote the L-function obtained by twisting Le{s) quadratically. 
Here d is a fundamental discriminant, i.e., d € Z — {1}, s.t. \ d for all odd primes 
p and d = 1 mod 4 or d = 8, 12 mod 16, and Xd is the Kronecker symbol. Then 
the twisted L-function (which is itself the L-function associated with another elliptic 
curve Ed) is given by 

n=l P f J 

where V'M is the principal Dirichlet character of modulus M: 

I U otherwise. 
The functional equation of this L-function is 
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In order to derive the 1-level density of the zeros near the critical point s = 1/2 
of L-functions in this family of quadratic twists, we consider the average over the 
family of a ratio of L-functions evaluated at different points: 

Xd{-M)uiE=+l 

This is an average over those twisted L-functions that have even functional equations 
and Q < d < X . Requiring an even functional equation imposes a restriction on d 
mod M. We follow the recipe of |CFK+05| . |CFZa| and the calculations in |CS07] to 
derive a formula for RE{ct, 7) via the ratios conjecture. Note that arriving at a ratios 
conjecture entails applying a list of manipulations, several of which introduce errors 
large enough to be significant. The miracle is that these errors appear to cancel 
out and the recipe yields formulae that have been checked numerically and against 
specific known cases in many different situations (see |CFZal ICSOTj ) . Recent work 
of Steven J. Miller [Milbj has shown that a rigorous calculation of the 1-level density 
for the family of real quadratic Dirichlet L-functions matches exactly, for a suitably 
chosen test function, the prediction obtained by applying the ratios recipe. See also 
[Sto] for further investigations of the ratios conjecture and the 1-level density of the 
same family of Dirichlet L-functions and [Milaj for Miller's extension of fMilb] to 
families of cuspidal newforms. 

We use (|2.3p to replace LE{s,Xd) in the denominator of (|2.6p by 



LE{s,Xd) ^ n 



where /U£;(n) is a multiplicative function defined as 

^iE{n) 



— A(p), if n = p 

V'Af(p),if n = p2 (2.8) 
0, if n = p^, k > 2. 

We use the approximate functional equation for the L-function in the numerator of 
LEil/2 + a,Xd) = ^ WTa)^ nV^- 

m<x \ / "<!/ 

+ remainder, (2-9) 

where M is the conductor of the elliptic curve E and xy = d"^/ {2-k). Therefore using 
the first sum of the approximate functional equation ()2.9p we get 

pi / ^ V- >^ X{m)nE{h)xd{mh) 

i?^(a,7):= 2^ 2^ ^i/2+a/,i/2+7 • (2.10) 

0<d<X h,m ^ 
Xd{-M)uOE=+^ 
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We denote by R^{a, 7) the expression that results from using the second sum in the 
approximate functional equation (j2.9p . Thus 

(2.11) 



Rsia, 7) ~ Rsia, 7) + ^|(a, 7)- 



The ratios recipe now cahs for a replacement of Xd{fnh) with its average over the 
family (the set of d's being summed over). We set 



0<d<X 
Xd{-M)ujE=+i 



1 and 



E 



1 



(2.12) 



0<d<X 
d=b mod M 



as the number of fundamental discriminants below X that we are summing over and 
note (see |CFK+05] . Theorem 3.1.1) 



- y 

X* „ 



Xdin) 



0<d<X 



Xb{g)0'{n) if n = g'D, with (DjM) = 1 and if all prime 
factors of g are prime factors of M 
otherwise, 



d=b mod M 

where 



an 



n 



p 



p + l 



(2.13) 
(2.14) 



This is to say that terms not of the form n = gO can be disregarded (this is the 
so-called 'harmonic detector' which is mentioned in jCFK"'"05] ). Since we are consid- 
ering only curves with prime conductor M, g is simply a power of M. Note that in 



M because d has been chosen such 
Xd{—^) since we are considering 



the cases we are interested in Xbid) = ^'^^ 9 
that Xb{M) = Xd{M) = iOE (we have Xd{M) 
only positive d). 

Concentrating on R^, we replace Xdi^h) with the average given by (j2.13p and 
so restrict the sum as follows: 



hm=DM<^ 



\{m) iJLE{h)a{'mh)uj\ 

mV2+a/il/2+7 



(2.15) 



with (□, M) = 1 and g divisible only by primes dividing M. We write this sum as an 
Euler product (for convenience denoting by h the exponent on primes dividing h in 
the sum above and similarly for m) and note that if m+h > 1 then a{p"^~^^) = p/p+1 
for primes not dividing the conductor, whereas a{p'^~^^) = 1 if the prime does divide 
the conductor. So we obtain 



i?lj(a, 7) «X*\^(a, 7)^1(0,7) 



(2.16) 
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where 

p\M \ ^ m,h>0 ^ / 

m+h>0 
m+h even 

^(«'7):=ii 2^ ^„.(i/2+a)+fe(i/2+-y) • (2.18) 

p|M \ /i,m>0 / 

Since ple{p^) = for most powers of p, we only need to consider /i = 0, 1, 2 in the 
sum in (j2.17p and /i = 0, 1 in (j2.18p . Then the Euler products become 



r) _|_ X \ ■'^-^ pm{l+2a) pl+Q+7 „m(l+2o) 
pfM \ \m=l ■'^ m=0 ■''^ 



1 Ad?-^™) \ \ 

'^pi+¥; „m(l+2a) 1 1 (2.19) 

m=0 



and 



^(«'7) = li ( 2^ 1 - pm(l/2+a)+l/2+y ) j" (2-20) 



p\M \m=0 



We now factor out the divergent part of using the Riemann zeta function and 
also, for convenience, we will factor out the symmetric square L-function associated 
with Le- This leaves us with a convergent Euler product. In the following, for 
simplicity, we shall only deal with elliptic curves with prime conductor, M. Recall 
that the Euler product of a Hasse-Weil L- function Le{s) coming from the elliptic 
curve E, with Dirichlet coefficients A(n) normalised so that the functional equation 
relates s to 1 — s, has the form 

Le{s) = - Kp)p-n-^ 11(1 - \{p)p-^ +p-^T'- (2.21) 

p\M p\M 
Now we can write this product as 

Le{s) = - «(p)p-^)-'(l - mP'T^ (2.22) 
p 

where 

a{p) + (5{p) = \{p) (2.23) 

and 

, , fo for p\M 

a{p)m = ' ^ 2.24 

I 1 for p f M. 
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Let LE(sym^,s) denote the symmetric square L- function. Then by definition (see 
|lwa97| . page 251) 

Ls(sym2, s) = - a\p)p-n-\l - a{p)mP'T'a " P\p)p~T'- (2-25) 
p 

We have (see |Con05j . page 236) 

A(m)A(n) = ^ \{mn/(f), (2.26) 

d\{m,n) 
(d,M)=l 

(where M is the conductor of E) and in particular we have for p \ M 

Kvf = A(/) + l (2.27) 
A(/™+i)A(p) = A(/'"+2)^^(p2m)^ (2.28) 

We wish to write the Euler product in (|2.25p in terms of A(p), so we start by using 
(f2:23]l to obtain 

2 , rrA Kp?-(^{p)m , «W(p)(A(p)'-«(p)/3(p)) 

Ls(sym,s) = [[[I + 



We now distinguish between p\M and p f M, and so, using (|2.27|) and (|2.24|) . we 
have 

MW.)^n(-^)"nO-f^-^-^;^)"'. (-0) 

p\M ^ -f^ / pI^j \ F F F / 

Now we reconsider the Euler products in (12.190 and (I2.20|) . In constructing 
ratios conjectures we usually allow —\< Rea < \ and logX <C Re7 < |, where 
the bounds at ^ allow us to control the convergence of Euler products of the type 
(j2.19p . In fact, in this application the real parts of a and 7 can be considered as 
very small. Thus we can write 

= n 1^1 I ^ is- ^(^'"^ V ^(^'"^'^ 

' ' ' \m=l m=0 



1 A(p^-) \ \ 

pm{l+2a) ) ) 
m=0 / / 



^/ A(p^) A(p^) + 1 1 . 

lll^+„l+2a + „l+27 + j' ^^■■^^> 
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where the • • • indicate terms that converge hke when a and 7 are small. We now 
use the following approximations to factor out the divergent or slowly converging 
terms. By ()2.30p we have 

Ls(sym2, 1 + 2a) = n (1 + ^ + • • • ) (2-32) 



and 



L£;(sym2,l + Q + 7)(;;(l + Q + 7) ^ 



Also, since there is only one prime that divides the conductor M, a factor of ^(1 + 27) 
will account for the divergence of the term ^i+i^ in (I2.3ip . 
Hence we can write 

V^ia,-f)V\{a,j) =YEia,-f)AE{a,j), (2.34) 

where 

V f ^ C(l + 27)^i;(sym^l + 2a) 

= C(l + a + 7)MsymM + a + 7)- ^''''^ 

^£;(a,7) is given by 



pfM \ \m=l ^ 

„m(l+2a) r,l+27 „m{l+2a) ) ) {^-M) 



pl+a+'y 



p|A/ \m=0 



A(p-)a;^ A(p) A(p™)a;™- 



p»n(l/2+a) pl/2+7 pm{l/2+Q) 

and is analytic as 0,7 ^ 0. Hence, by recalling (|2.1U|) . we find 



^^(«,7)« Yl YE{a,^)AE{a,-f). (2.37) 



0<d<X 
Xd{-M)uJE=+l 



We obtain the other sum B?^{a, 7) in ()2.1ip by using the second term in the approx- 
imate functional equation (j2.9p and carrying out exactly the same steps as above: 

^l(«'7)- E ^^^r^yE{-an)AE{-an)- (2.38) 



27r / r(l + a 
By applying the ratios conjecture recipe, we therefore have the result: 



Q<d<X 
Xd{~M)uiE=+^ 
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Conjecture 2.1 (Ratios Conjecture). For some reasonable conditions such as —\ < 
Rea < J, i^Jjf ^ 1^67 < j and Ima,Im7 <^ X^~^, we have 

0<d<X ' 
Xd{-M)ujE=+l 

Xd{-M)uiE=+l 

where Ye and Ae are defined at \2. 35\) and ^2.36\) . respectively, M is the (prime) 
conductor of the L-function Le{s) and loe is the sign from its functional equation. 

We note that the error term 0(X^/^+^) is part of the statement of the ratios 
conjecture; the power on X is not suggested by any of the steps used in arriving at 
the main expression in Conjecture 12.11 At the end of Section [3] we propose that the 
hmited data we have available supports a power saving on the error term, but not 
necessarily a power of 1/2. 

To calculate the 1-level density we actually need the average of the logarithmic 
derivative of L-functions in this family, so we note that 

oh^x LE{l/2 + r,Xd) da ' ^ > 

Xd{-M)wE=+l 



a=-y=r 



Using ()2.28p for primes not dividing M and the multiplicativity of \{p) for p\M, we 
get AE{r,r) = 1 and we have, with 

AUr,r) = ^AE{a,j)\^^^^^, (2.40) 

d , , ^ C'(l + 2r) , , , L'p(sym2,l + 2r) , , , .i . x 

^YEAEia,,)l^^^^ = - '-^^AE{r,r) + ^ J^ AE[r,r) + A\ir,r) 

_ C{l + 2r) L-^(symM + 2r) 

" C(l + 2r) + LE(sym2, 1 + 2r) + ^^'^^^ 



da 



and 




2a 

^^i^{y^(-a,7)^E(-a,7)} 

1 (1 + aj a=7=r 

r(l - a) C(l + 2r)L£;(sym2, 1 - 2r) 



r(l + a) Ls(sym2,l) 
Therefore we have for the logarithmic derivative the following: 



AEi-r,r). (2.42) 
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Theorem 2.2. Assuming the Ratios Conjecture \2.1\ and i^gX ^ R'e(r) < ^ and 
Im(r) <C X^~^ , the average of the logarithmic derivative over a family of quadratic 
twists (with even functional equation) of the L-function of an elliptic curve with 
prime conductor M is 



L'j,{l/2 + r,Xd) 

Xd(-M)ujE=+l 



^ LE{ll2 + r,Xd) 

0<d<X ' 




[ C'(l + 2r) L' (sym2,l + 2r) .^ , , , , 

Xd{-M)LOE=+l 

r(l-r)C(l + 2r)L^(symM-2r) \ 
r(l + r) L^(sym2,l) r,r)j+0[X ). 

HereuiE is the sign from the functional equation of Le, LE{sym.'^ , s) is the associated 
symmetric square L-function (defined at i2.25\) ). and Ae and A\, are arithmetic 
factors defined at ^2.36\) and i2.40\ ), respectively. 

Let 7d denote the ordinate of a generic zero of LE{s,Xd) on the half line. We 
consider the 1-level density 

^i(/)^= E (2.44) 

0<d<X 7d 
Xd(-M)iJE=+l 

where / is some nice test function, say an even Schwartz function. By the argument 
principle we have 

^^(/)- E ^([-L .)|^/H(-V2))^« (2.45) 



Xd(-M)a;B=+l 

where (c) denotes a vertical line from c — zoo to c + ioo and 3/4 > c > 1/2 + 1/ log X. 
The integral on the c-line is 

2vry_oo L£;(l/2 + (c- 1/2 + zt),xd) 

Xd(-MVB=+1 

The sum over d can be replaced by Theorem 12.21 The bounds on the size t coming 
from the ratios conjecture should not limit us here. It is not entirely known in 
what range of the parameters the ratios conjecture holds, but the test function / 
can be chosen to decay sufficiently fast that the tails of the integrand, where the 
ratios conjecture might fail, will not contribute significantly. (See the 1-level density 
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section of |CS07j for more detailed analysis.) Next we move the path of integration 
to c = 1/2 as the integrand is regular at t = and get 

2t^J-J^' \ C(l + 2it) Li5(sym2,l + 2it) . > \ ) 

Xd{-M)uJE=+l 

■ Midi V" r(l - it) C(l + 2it)L.(.ym^ 1 - ^A^i-^t.^Xt 



\ 2tt J T{l + it) Lij(sym2,l) 

For the integral on the line with real part 1 — c, we use the functional equation 
LEis, Xd) = Xd[-M)uEX{s, Xd)LE{l - s, Xd) (2.48) 



with 



to obtain 



L'eC^ - s, Xd) _ X'{s, Xd) L'j^{s, Xd) 



LE{^-s,Xd) X{s,Xd) LE{s,Xd)' 
The logarithmic derivative of ()2.49p evaluated at s = 1/2 + a is 



(2.50) 



X' {1/2 + a,xd) ( ^\d\ \ T' V 

For the integral on the (1 — c) line we change variables s — > 1 — s and use ()2.50p . 
We thus obtain finally the following: 

Theorem 2.3. Assuming the Ratios Coniecture \2.1\. the 1-level density for the zeros 
of the family of even quadratic twists of an elliptic curve L- function Le{s) with prime 
conductor M is given by 

0<d<X 7d 
Xd{-M)LUE=+l 

^ f°° . V- ( ( VM\d\\ V' , T\ 

0<d<X \ \ / 

Xd{-M)ujE=+l 

,2 



+ 2 



C{l + 2it) LE{sym^,l + 2it) e\ , ) V ; 

dt 



M\d\\ r(l-it)C(l + 2it)L£;(sym2,l-2it) 



I 27r / r(l + ii) Lij(sym2,l) 
+ 0(Xi/2+^), 
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where jd is a generic zero of Le{s, Xd), f is an even test function as described above, 
oje is the sign from the functional equation of Le, LE{sym^ , s) is the associated 
symmetric square L-function (defined at ^2.25\) ). and Ae and are arithmetic 
factors defined at i2.36\) and i2.40\ ), respectively. 

3 Numerical test 

We test our prediction - namely formula ()2.52p - for the 1-level density with a con- 
crete example (see figure [2]). We pick the elliptic curve En with (ai, 02, 03, 04, ae) = 
(0,-1,1,0,0) in the WeierstraB form 

+ aixy + a^y = + 02^^ + 04X + (3.1) 

giving 

En:y^ + y = x^- (3.2) 

and consider the even quadratic twists of its associated L-function with fundamen- 
tal discriminants between and 40,000. We are interested in the 1-level density of 
unsealed zeros from up to height 30. The numerical data is obtained from Ru- 
binstein's Icalc [Rub]. In the range considered we find 11,135 quadratic twists, of 
which 5,562 are even ones with a total of about 590,170 zeros. In figure[2]we obtain 
the solid curve from the histogram of this zero data by choosing a binsize of 0.1 
and dividing by both the number of quadratic twists with even functional equation, 
and the mean density of zeros log(\/ITX/(27r)). 593 of the L-functions with even 
functional equation have (at least) a double zero at the central point; these zeros 
at the central point are not plotted in figure [2 The dashed curve is obtained from 
the formula ([TS^ with X = 40, 000 and f{t) = 5{t - x) + 6{t + x) for x between 
and 30. This curve is scaled like the data curve by dividing through by the number 
of quadratic twists with even functional equation and the mean density of zeros. It 
was computed using a combination of Mathematica and C++. The coefficients \{p) 
appearing in the arithmetic factor AE{a, 7) were computed using PARI. To compute 
coefficients of prime powers A(p™) for p \ M the following recursion formulas (see 
|HM070 were used 

^ Mpf' - g - ,)) A(p-) (3.3) 

A(p^™«) = A(p)^"'«-g(fr-r)" (J- X'l)) <^-'" 

In general there is good agreement between the data and the theoretical curve, 
which captures the main features of the data. We would expect better agreement 
with a larger set of data, since the data seems not yet to have resolved all the peaks 
further out along the axis. 
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Figure 2: 1-level density of unsealed zeros from up to height 30 of even quadratic 
twists of Lej^i with < d < 40, 000: prediction (dashed), from (|2.52|) . versus numer- 
ical data (solid) 
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A closer look reveals that the 1-level-density is strongly governed by the non- 
trivial zeros of C(s) and L(sym^,s): we observe that some dips of the data curve 
are located at 7/2 where 7 is the ordinate of a non-trivial zero of the Riemann zeta 
function. This is captured in the term 

- 4^ (3.5) 

C(l + 2it) ^ ' 

of our conjecture for Si{f). In figure [3] we mark the position of a non-trivial zero of 
the Riemann zeta function on our conjectural answer by a *. These * are all localised 
in or around a neighbourhood of a dip. This phenomenon has been encountered 
before, in the study of lower order terms of the number variance [Ber88] and the 
correlation functions |BK99[ lBK96al ICSOTl [CM ICS] of the Riemann zeros, and in 
the one- level density of other families of L-functions |CS07j . 

On the other hand we observe that some peaks are located at 7/2 where 7 is the 
ordinate of a non-trivial zero of L£;(sym^, s). This is captured in the term 

L'j,{sym\l + 2it) 

LE{sym'^,l + 2it) ^ ' ' 

of our conjecture for Si{f). In figure [3] we mark the position of a non-trivial zero 
of L£;(sym^, s) by a o. The majority of these os are localized in or around a neigh- 
bourhood of a peak. In particular, we observe that if a zero of the Riemann zeta 
function is close to a zero of L(sym^,s) then these zeros are localised in or around 
a dip. Hence, zeros of the Riemann zeta function seem to dominate the behaviour 
of the 1-level-density more than the zeros of -L(sym^,s). This may be explained 
because the density of the Riemann zeros in this range is smaller than that of the 
zeros of L(sym^, s) and so in terms of the mean zero density the one-line is closer to 
the half-line in the case of the Riemann zeta function. Therefore one would expect 
the Riemann zeros to have a larger effect. 
The term 

M\d\\ T{l-it)C{l + 2it)LE{symM-2tt) 

AE[-it,it), (3.7) 



2tt J T{l + it) Li5(sym2,l) 



from ()2.52p . makes its most obvious contribution by causing the oscillation near 
the origin of the plot of our conjectural answer for the 1-level density. The factor 

results in oscillations on the scale of the mean density of the zeros of 



2tt 

the original L-function, Le- 

In summary, we notice that the lower order terms dominate the behaviour of 
the zeros when we are far from the limit of infinite conductor (in the family of 
quadratic twists, E^, the conductor increases with d). This becomes more obvious 
when we compare our conjectural answer for finite conductors with the limiting 
theoretical result: in figured] we consider the scaled 1-level density of S0(2N) in 
the limit — > 00 against our conjectural answer (also scaled) for finite conductor. 
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Figure 3: Effects of non-trivial zeros of the Riemann zeta function (indicated by *) 
and the non-trivial zeros of -L(sym^,s) function (indicated by o) on the conjectural 
formula ()2.52p for the 1-level density of unsealed zeros from up to height 30 of 
even quadratic twists of Lei^ with < d < 40, 000 
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We observe convergence to the limiting theoretical result as we increase X, the cut- 
off point for d. The observed effects of the arithmetical terms for small and finite 
conductors are washed out and shifted away from the origin in the large conductor 
limit. 

To further understand the approach to the limiting distribution, we calculate the 
1-level density for scaled zeros and recover the limit and the next to leading order 
term from (|2.52|) . As a first step we rescale the variable t in (|2.52|) as 



and define 
where 



L := lo; 

and get, after a change of variables. 



r = t(L/7r) 

fit) = git{LM), 
MX\ 



27r 



(3.8) 
(3.9) 
(3.10) 



0<d<X 7d 
Xdi-M)u;E=+l 



1 

2L 



E 



0<d<X 
Xdi-M)u}E=+l 



2 log 



27r 



r' / iiTT 



r' / iTTT\ 



+ 2 



2iTTT \ 



+ 



C(l + ^) ■ LE{sym\l + ^ 
M\d\ \ -2W^ r(l - Y) C(l + ^)Lg(sym2, 1 - ^ 



\ L ' L J 



2tt / 



r(i + Y) 



L£;(sym2, 1 

, / ivrr ivrrx 

"^4-— '— ) 



(3.11) 



+ 0(Xl/2+=). 

We write the number of fundamental discriminants less than or equal to X as 

X* := 1. (3.12) 

0<d<X 
Xd{~M)uJE=+l 

Using the Euler-Maclaurin formula we make the approximation 

\/mx\ 



E 



log 



0<d<X 
Xd{-M)uJE=+l 



. 2- i 



X* 



log 



27r i 



+ 0(x^/^+'). (3.13) 
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Figure 4: Scaled limiting 1-level density of S0{2N) (solid) versus scaled formula 
(fXTT]) divided by X* (dashed) for: < d < 40,000 (top left), < d < 10^ (top 
right), < d < 10^° (middle left), < d < 10^° (middle right), < d < 10^° 
(bottom left), < d < 10^°° (bottom right) 
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In the same manner we have 



/ \ —2iTTT/L 

J2 =X* (1 + ^ + 0(L-2))e-2-- + 0(XV2). 



0<d<X 
Xd{-M)^E=+l 



2tt ] ^ L 

(3.14) 



Writmg 



oo 



1 , , ( _1 1" 

C(s + 1) = - + J;^^7ns^ (3.15) 

we have 



s ^ n\ 

rt=0 



7^ = - + 7 + (-7^ - 27i)« + 0(.2), 
C(l + s) 



(3.16) 



where 7 = 70 is Euler's constant, and so 



and 



C(i + gf^) ^ L 



+ j + 0{L-'). (3.18) 



Simple Taylor expansions of the other factors in ()3.1ip lead us to, with the relation 
between / and g given in (13. 9p . 



X* X* ^ vr 

Q<d<X -Id 



where 



and 



with 



/■°° , sin(27rT) 1 + cos(27rT) 7rTsin(27rr) ^/l\\, 



a.^l + 27-^kO,0)-||^ (3.20) 



«2 = 2 + 47 + 37' - 271 + B'{0) + 2^B'{0) - 2 ^/^'^'^l']\ 

L(l) Li=;(sym2,l) ^ 4 ^Ls(sym2,l)' ^ ' > 



B'{0) = ^AE{-r,r) and i?"(0) = -^^^(-r, r) . (3.22) 
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In order to obtain ()3.20p we use the following identity 




1(0,0). 



(3.23) 



We establish identity (j3.23p by simple algebra and using (|2.28p for primes not di- 
viding M, the multiplicativity of X{p) for p\M and AE{r,r) = 1. 

This work was initially conceived to investigate the unexpected numerical results 
found by Steven J. Miller jMilOG] near the origin of the histogram of the distribu- 
tion of the first zero above the central point of a family of rank zero L-functions. 
He observed very few examples of zeros lying close to the central point. That is, 
he observed the phenomenon of repulsion of zeros from the central point which we 
know, from rigorous work on the 1-level and 2-level densities |You06l IMil02l IMil04j 
does not persist in the large conductor limit. Since the 1-level density (a histogram 
of all zeros) and the distribution of the lowest zero (a histogram of the lowest zero 
of each L-function) are the same for very small distances from the central point, it 
is natural to enquire whether the ratios conjecture yields a formula for the 1-level 
density which would display and explain Miller's observed repulsion at finite con- 
ductor. Although it can be seen from figured] that the formula ()3.1ip is significantly 
smaller near the origin than the limiting curve, and approaches it from below as the 
conductor increases, there is no evidence of repulsion. This is a major discrepancy 
from the data, as seen in figured) away from the critical point we have a nice match 
between the prediction and the data while near the critical point we find fewer zeros 
in the data than predicted by our formula. It is most interesting that the main 
terms of the ratios conjecture do not capture this important feature. Of course, the 
natural question is whether this contradicts the ratios conjecture, or whether the 
discrepancy can by accounted for by the error term. As expected due to the limited 
data available, the test described below is inconclusive, but shows signs that the 
error term in the ratios conjecture (and hence on the one level density in (|2.52p ) is 
of the form X^"*"^, for b < 1. The ratios conjecture is usually stated with b = 1/2. 

We fix several sample points at various distances away from the critical point 
and measure the difference between the main terms of our prediction (that is, the 
sum over d inside the integral in (j2.52p ) and the data. In fact, we compare the 
normalised versions of our prediction and data by dividing through by the number 
of fundamental discriminants X* less than X and the mean density of zeros. So let 
us denote this difference between the main terms of the normalised theory and the 
data at a fixed height t and fixed X by A(t, X). Since we have divided by X* , which 
is proportional to X, this difference is expected to be of size 




(3.24) 



The quantity we will plot is 




(3.25) 
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and if the ratios conjecture with error term X^^'^ is correct then we would expect 



=6-1 + 0(^^1^) (3.26) 

as X — > oo. 

In figure [5] we plot the quantity QA{t-,X) for < X < 400,000 and for various 
fixed sample points = 0.01, t2 = 0.02, is = 0.03, t4 = 0.04, ts = 0.05, te = 0.4 
and, tj = 0.6. We notice that the curves are much smoother for sample points 
near the critical point, i = 0, eg. ti,t2,t3. In the range < X < 400,000 these 
points are well inside the region where the zero data shows repulsion at the critical 
point; see figured! Thus the difference between the theory (smooth curve in figure 
[T]) and data (histogram) does not change sign as X increases. Presumably it is the 
amplification of such sign changes by the logarithm in (j3.25p that is responsible for 
the jagged curves in figure [5] for sample points t^jt^ and tQ. 

We see also that the curves at sample points close to the critical point appear 
at first sight to indicate a larger error term - in fact, over this range of X the ti 
curve implies 6 — 1 > 0! If a limit such as (j3.26p exists, it does not seem to behave 
uniformly in t. However, the ti,t2 and curves are decaying as X increases and we 
do not have enough data to see what their final behaviour will be. We remember 
that the convergence is like log log X/ log X, so we would need much more data to 
be able to make a sensible conclusion about the size of the error term. 

Also, it is interesting to note that at the right hand side of figured the = 0.03 
curve has decayed to a level comparable to the curves of the sample points that are 
more distant from t = 0. Examining figure [H it appears that the area of major 
discrepancy between the ratios conjecture prediction and the data (that is, where 
the data shows repulsion from the critical point at t = 0) lies between t=0 and about 
t = 0.03. We expect that this region will narrow as the range of discriminants, d, 
increases, and this is born out by comparing the two pictures in figure [U the data 
grows more quickly to the height of the solid curve in the right hand picture where 
< d < 400, 000, than in the left hand picture where < d < 100, 000. Thus at 
the right hand edge of figure O the point ^3 = 0.03 is about to move into the region 
where there is good agreement between the ratios conjecture prediction and the 
data. Making a speculative conclusion from the limited data available, this suggests 
that the curves for ti and t2, or any other fixed t, would also decay to this level if we 
could gather enough data to shrink the area of discrepancy at the origin of figure [J 
to a narrow enough band. 

It is impossible to say from the available data what the exponent b in the error 
term of the ratios conjecture is. There is certainly no evidence to suggest b = 0.5, 
but the possibility that the curves in figure [5] would decay to —0.5 if we could vastly 
extend the rage of the plot is not ruled out. However, figure [5] certainly appears to 
suggest that 5 < and so the error term is a power of X smaller than the main 
term. 
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Figure 5: discrepancy QAit^X) from ()3.25p between prediction and data for the 
1-level density of even quadratic twists of with < d < 400, 000. 



4 Summary 

We find that the ratios conjecture provides a formula for the one level density of 
zeros of a family of quadratic twists of an elliptic curve L-function that agrees with 
data for finite conductor, except in the vicinity of the critical point, t = 0, and 
explains the arithmetic nature of the lower order terms which entirely dominate the 
behaviour of the statistic away from t = 0. The ratios conjecture prediction, when 
properly scaled, approaches the limiting SO{2N) random matrix result as the family 
of elliptic curves includes those with larger and larger conductor. This supports all 
the available evidence that SO{2N) is the correct limit for zero statistics in this 
family. It is very interesting that the ratios conjecture prediction does not capture 
the phenomenon of zero repulsion from the critical point, t = 0, but the data we 
have available certainly allows for the ratios conjecture to be correct with some 
power 6 < 1 of X in the error term; the discrepancy between the ratios conjecture 
prediction and the data (at the origin of figured]) can quite possibly be contained 
in the error term. 

In ongoing work of the authors in collaboration with E. Dueiiez and S. J. Miller 
we propose an explanation for the observed repulsion of zeros near the central point 
for finite conductor and a random matrix model that captures the phenomenon. 
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